An efficient meshfree point collocation moving least squares method to solve the interface problems with nonhomogeneous jump conditions

We are going to study a simple and effective method for the numerical solution of the closed interface boundary value problem with both discontinuities in the solution and its derivatives. It uses a strong‐form meshfree method based on the moving least squares (MLS) approximation. In this method, for the solution of elliptic equation, the second‐order derivatives of the shape functions are needed in constructing the global stiffness matrix. It is well‐known that the calculation of full derivatives of the MLS approximation, especially in high dimensions, is quite costly. In the current work, we apply the diffuse derivatives using an efficient technique. In this technique, we calculate the higher‐order derivatives using the approximation of lower‐order derivatives, instead of calculating directly derivatives. This technique can improve the accuracy of meshfree point collocation method for interface problems with nonhomogeneous jump conditions and can efficiently estimate diffuse derivatives of second‐ and higher‐orders using only linear basis functions. To introduce the appropriate discontinuous shape functions in the vicinity of interface, we choose the visibility criterion method that modifies the support of weight function in MLS approximation and leads to an efficient computational procedure for the solution of closed interface problems. The proposed method is applied for elliptic and biharmonic interface problems. For the biharmonic equation, we use a mixed scheme, which replaces this equation by a coupled elliptic system. Also the application of the present method to elasticity equation with discontinuities in the coefficients across a closed interface has been provided. Representative numerical examples demonstrate the accuracy and robustness of the proposed methodology for the closed interface problems. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1031–1053, 2015     

Adaptive discretization for signal detection in statistical inverse problems

We discuss statistical tests in inverse problems when the original equation is replaced by a discretized one, i.e. a linear system of equations. Previous studies revealed that using the discretization level as regularizing procedure is possible, but its application is limited unless discretization is restricted to the singular value decomposition, see C. Marteau and P. Mathé, General regularization schemes for signal detection in inverse problems, 2013. General linear regularization may circumvent this, and we propose a regularization of the discretized equations. The discretization level may be chosen adaptively, which may save computational budget. This results in tests which are known to yield the optimal separation rate up to some constant in many cases.     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.     

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.     

Du Fort-Frankel格式及DFF-I并行格式的稳定性

 本文对抛物型方程的Du Fort-Frankel(DFF)格式以及基于该格式构造的并行差分格式(DFF-I)进行了稳定性分析。采用矩阵分析方法, 证明了其无条件(LR)稳定性, 给出了DFF格式的稳定性系数的最小值的上界估计, 结果表明其与网格比有关, 从而DFF格式并非绝对稳定。本文改进了并行差分格式(DFF-I)的稳定性分析结果, 证明了其增长矩阵的谱半径严格小于1, 从而具有长时间稳定性。数值算例验证了DFF-I格式具有空间二阶精度, 且有很好的稳定性。     

液滴撞击加热壁面传热实验研究

本文采用高速摄像仪对水滴和乙醇液滴撞击加热壁面后的蒸发过程进行了实验观测, 分析了液滴撞击加热壁面后的蒸发特性参数. 实验中, 两种液体初始温度均为20 ℃, 不锈钢壁面初始温度范围为68-126℃. 水滴初始直径为2.07 mm, 撞击壁面时Weber 数为2-44; 乙醇液滴初始直径为1.64 mm, Weber数为3-88. 结果表明, 液滴受到重力、表面张力及流动性的影响, 在蒸发过程的大部分时间内, 水滴高度持续降低而接触直径几乎不变; 蒸发后期, 液滴发生回缩, 水滴的接触直径、高度和接触角出现振荡现象. 乙醇液滴的接触角随时间的增加呈现先减小随后保持不变的趋势, 而接触直径和高度则持续减小, 直到液滴完全蒸发. 液滴蒸发总时长与液体物性和壁面温度有关, 随壁面温度的升高而减小, 与液滴撞击壁面时的Weber 数无关. 同时, 随着壁面温度的升高, 液滴显热部分占总换热量的比重增大, 显热部分能量不可忽略, 本文实验条件下得到水滴的平均热流密度为0.014-0.110 W·mm^-2.     

立体发光LED灯片的光学性能研究

立体发光灯片的关键是在保证光通量的前提下,如何获得均匀的正反面发光.本文对3种荧光粉平面涂覆立体灯片的正面和反面的瞬态和稳态光通量、色温、色品坐标进行了研究,发现直接涂覆荧光粉胶的灯片总的光通量最大,但正反面的发光性能相差很大;在涂覆荧光粉之前先平面涂覆一层硅胶,可获得正反面发光较均匀的灯片;而在该硅胶层内掺杂扩散粉的灯片,不但没有提高灯片的正反面发光均匀度,反而大幅度降低了光通量值.此外,测试3种样品的发光角度,发现其发光角度分布相似,而发光强度与积分球测试的光通量比值相近.结果表明,涂覆荧光粉之前先平面涂覆一层硅胶工艺,既可以保证立体发光灯片具有较大的光通量又可以获得均匀的正反面发光.     

��ͨ����Һ̬����ٶ���������Һ̬��¶ȷֲ���Ӱ��

应用商业软件ANSYS CFX计算了等离子体热通量和液态锂流速对自由流动液态锂温度分布的影响。计算结果表明,导向槽中心附近液态锂温度较高,冷却水入口和出口对应位置液态锂温度最低。液态锂出口温度随着等离子体热通量的增大而线性升高,冷却水流速为1.5m·s-1,热通量分别为0.1MW·m-2和1MW·m-2时,液态锂在出口处对应的温度分别为255.3°C和458.6°C。增大液态锂流速,导向槽内液态锂的温度逐渐降低,但温度变化的幅度较小。计算结果对液态锂回路安全稳定运行提供了一定参考。     

Flux calculation at the interface between two rock types for two-phase flow in porous media

When simulating two-phase flow in porous media, one has to consider the case where there is a discontinuity in the medium. There relative permeabilities and capillary pressure functions may change and we address the problem of calculating the convective part of the numerical flux at the interface between the two rock types. Several solutions are compared.     

一种容积法测小流量研究

针对高压汽水两相流实验中小液量测量困难这一问题，运用重力压差主要由液相密度决定的原理，提出了由圆管内重力压差变化线性分求流量的方法，实验表明这一方法线性度很好，本文对这一方法进行了理论分析，得到了计算修正公式，同时评估了计算带来的误差。